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Abstract —Application of filter bank multicarrier (FBMC) as 
an effective method for signaling over massive MIMO channels 
has been recently proposed. This paper further expands the 
application of FBMC to massive MIMO by applying frequency 
spreading equalization (FSE) to these channels. FSE allows us 
to achieve a more accurate equalization. Hence, higher number 
of bits per symbol can be transmitted and the bandwidth of 
each subcarrier can be widened. Widening the bandwidth of 
each subcarrier leads to (i) higher bandwidth efficiency; (ii) 
lower complexity; (lii) lower sensitivity to carrier frequency offset 
(CEO); (iv) reduced peak-to-average power ratio (PAPR); and 
(Iv) reduced latency. All these appealing advantages have a direct 
impact on the digital as well as analog circuitry that is needed 
for the system implementation. In this paper, we develop the 
mathematical formulation of the minimum mean square error 
(MMSE) FSE for massive MIMO systems. This analysis guides 
us to decide on the number of subcarriers that will be sufficient 
for practical channel models. 

I. Introduction 

Based on the recent discussions on the 5**^ generation of 
wireless communication networks (5G), massive MIMO, as 
a candidate, can play a pivotal role in these networks, HI. 
Additionally, as a fundamental building block in the physical 
layer of 5G networks, the need for signaling techniques that 
are compatible with the emerging applications such as the 
machine to machine (M2M) communications and Internet of 
Things (loT) is greatly emphasized ||2l. Filter bank multicarrier 
(FBMC), a method that was initially proposed about 50 years 
ago, El, El, is recently being considered as a candidate for 
5G systems, Cl. 

Massive MIMO is a code division multiple access (CDMA) 
like system that may be used to increase the capacity of 
multiuser networks, 0. As it is shown in 0, by increasing 
the number of antennas at the base station (BS), the effects 
of noise and multiuser interference (MUI) start to vanish until 
they will be completely removed as the number of BS antennas 
tends to infinity. In the recent works 0 and 0, FBMC was 
proven to be a good match for massive MIMO and vice versa 
as they can both bring a great number of appealing properties 
into the picture of 5G systems. From network capacity view¬ 
point, this combination is of a paramount importance as the 
same spectrum is not only being reused by all the users (an 
advantage of massive MIMO) but it is being used in a more 
bandwidth efficient way. 

A few recent publications have introduced an alternative to 
the traditional polyphase-based implementation of the FBMC 
systems that is termed as frequency spreading FBMC (FS- 


FBMC), 0- lfT3ll . This technique offers a number of advan¬ 
tages that put FBMC in a strong position as a candidate for the 
future wireless standards. FS-FBMC was explained for the first 
time in a document that came out of the PHYDYAS project, 
El. Further details of FS-FBMC were presented, by the same 
author, in ifTOI . The sensitivity of FS-FBMC to the timing 
offset is studied in ifTTIl where the robustness of the FS-FBMC 
frequency domain equalizer to the timing errors compared with 
its counterpart with a polyphase-based receiver structure was 
reported. Berg et ak, na, have explored the utilization of FS- 
FBMC to access TV white spaces. All these works have used 
the filter designed by Martin IH and Mirabbasi and Martin 
ifTSll for the implementation of FS-FBMC system. 

As an extension to our previous work, 0, in this paper, we 
introduce the frequency spreading equalization (FSE) concept 
to the realm of FBMC-based massive MIMO systems. Thanks 
to the better equalization capability of FS-FBMC systems, we 
will show that FSE brings a great amount of performance 
enhancement compared with the results of 0 . This allows the 
usage of even wider subcarriers than what was suggested in 0 
which in turn further reduces the latency due to the synthesis 
and analysis filter banks. Furthermore, this unlocks utilization 
of higher order modulation schemes in FBMC-based massive 
MIMO systems. Using a smaller number of subcarriers leads 
to a lower complexity, lower sensitivity to CFO, reduced PAPR 
and reduced latency. These advantages simplify the digital and 
analog circuitry that is needed for the system implementation. 

Widening the bandwidth of each subcarrier is equivalent to 
shortening the length of the prototype filter. This, in turn, has 
the impact of increasing the bandwidth efficiency of FBMC, 
as the ramp-up and ramp-down of the FBMC signal will be 
shortened. A factor that limits utilization of FBMC systems 
in some applications like bursty M2M communications is the 
long ramp-up and ramp-down of their signals. To tackle this 
problem, some researchers have introduced the concept of 
FBMC systems with circular pulse-shaping, m-iisi, where 
the prototype filter transients can be removed through the 
so called tail biting property. Even though these systems 
can remove the prototype filter transients from their signals, 
they need to add a cyclic extension to form their signal. 
Additionally, they suffer from higher out of band emissions 
than the FBMC systems with linear pulse-shaping, insi, ca. 
The approach in this paper is to reduce the ramp-up and ramp- 
down of the FBMC signal through widening the subcarriers 
while achieving a good equalization performance. This re- 


duced transient duration of the linearly pulse-shaped FBMC 
signals is equivalent to the cyclic extensions being used to 
form the circularly pulse-shaped FBMC signals. Therefore, in 
terms of bandwidth efficiency, both systems seem to be about 
the same while linearly pulse-shaped FBMC systems keep the 
advantage of having a lower out of band emissions. 

The rest of the paper is organized as follows. FBMC signal 
construction and the frequency spreading method are explained 
in Section |II] The proposed frequency spreading equalization 
for the massive MIMO systems is discussed in Section m 
where we have also derived analytical signal to interference 
plus noise (SINR) expressions. Various aspects of our design 
are analyzed through computer simulations in Section |IV] 
Finally, the paper is concluded in Section |Vl 

Notations: Scalars, vectors and matrices are represented in 
regular letters, boldface lower and upper case letters, respec¬ 
tively. In denotes an x TV identity matrix. The matrix 
or vector superscripts (■)^ and (•)^ indicate transpose and 
conjugate transpose, respectively. || • || and | • | represent Eu¬ 
clidean norm and absolute value, respectively. 5R{-}, S5{-}, E[-] 
and identify the real and imaginary parts of a complex 
number, expectation and inverse of a matrix, respectively. 
Linear convolution is denoted by ★, and the superscript * 
represents the complex conjugate. 

II. CMT AND ITS Frequency Spreading 
Implementation 

In cosine modulated multitone (CMT), which is a type of 
FBMC similar to the widely known OQAM-FBMC technique 
ED, real-valued data symbols are spread across time and 
frequency such that if the symbol time spacing, i.e., symbol 
period, is denoted by T, the symbol frequency spacing is 
F = Moreover, subcarrier phases are toggled between 0 
and tt /2 so that adjacent subcarriers have a quadrature phase 
difference. This phase-adjustment is necessary for the CMT 
receiver to be able to recover the transmitted data symbols 
free of interference along time and frequency. The hlter p{t), 
which is used to pulse shape the data symbols, is a square- 
root Nyquist hlter, i.e., q{t) = p{t)-kp*{—t) is a Nyquist hlter, 
where the spacing between its zero-crossings is equal to 2T. 

In practice, the CMT signal is efficiently synthesized at the 
transmitter and analyzed at the receiver using the polyphase 
network (PPN) ll22l . 12^ . However, as it has been recently 
suggested by Bellanger a, the frequency spreading structure 
can be utilized to gain a number of advantages compared 
with the traditional PPN method. In particular, an interest¬ 
ing property of FS-FBMC structure is its high-performance 
equalization capability. This equalization, namely, FSE, can 
be utilized in the context of massive MIMO as suggested 
in Section |III] In the rest of this section, and for our later 
reference, we review the frequency spreading structure when 
applied to CMT, which we hereafter call FS-CMT. 

The basic idea of the frequency spreading method is that a 
discrete-time square-root Nyquist filter, p[n], of size N = KL, 
can be synthesized in the frequency domain using a set of 


L/2 samples 

< -► 



Fig. 1. The overlap and add operation. 

2K — 1 distinct tones, as 

p(^) = X! (1) 

Here, Ck is a real-valued coefficient that indicates the fre¬ 
quency response of the p{n) at the frequency ojk = K is 
the overlapping factor, i.e., it signifies the number of symbols 
that overlap in the time domain and L is the sample spacing of 
the zero-crossings of the Nyquist filter q{n) = p{n)-kp*{—n). 

When the above idea is extended to the multicarrier case, 
the symbols at each subcarrier should be pulse-shaped in the 
frequency domain using the set of Cfc’s. An inverse discrete 
Eourier transform (IDET) block of size N can then be used 
to calculate the time domain samples from their frequency 
domain counterparts. Casting the above discussion into the 
mathematical form, the ES-CMT symbol can be obtained as 

x(n) = F]^A€>s(n). (2) 

Here, s(n) = [so(n) si(n) SL_i(n)]"^, where Sfc(n) 

indicates the data symbol at the subcarrier and n**' 
symbol time, indicates the phase-adjustment matrix which 
is an E X L diagonal matrix with the diagonal elements 
,... A is the spreading matrix of size 

N X L whose fc*** column contains the coefficients c^’s that 
are centered at the center frequency of the subcarrier, 
indicates the IDET matrix and x(n) is an A^ x 1 vector 
indicating the ES-CMT symbol in time domain. 

Noting that the size of each ES-CMT symbol is N and 
the symbol spacing is equal to L/2, the transmit signal can 
be constructed by performing the overlap-and-add operations 
on the symbol vectors as illustrated in Eig. [D Assuming an 
ideal channel and having a perfect synchronization between 
transmitter and receiver, the N samples that correspond to the 
symbol are collected in the vector f(n), at the receiver 
side, and the data symbols are recovered according to 

s(n) = 5R {$"^A'^FArf(n)} . (3) 

Since we assumed an ideal channel, equalization was not con¬ 
sidered here. In the presence of the channel, the equalization 
can be embedded after obtaining the DET output samples, m. 
In the next section, we extend the equalization scheme that was 
originally proposed for single antenna systems in ITOl to the 
context of massive MIMO. 
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III. MMSE-FSE IN Massive MIMO 

An interesting property of EBMC in massive MIMO chan¬ 
nels that is called self-equalization was recently reported in 
Q. Due to that property, a channel flattening effect takes 
place as a result of linear combining of the channel gains in 
different receiver antennas during the equalization procedure. 
As a result, wider subcarriers in frequency can be utilized 
which brings a number of benefits to EBMC systems in 
massive MIMO application. The benefits of EBMC-based 
massive MIMO systems are well articulated in Q and m. 
As it is shown in im, a better channel equalization can 
be achieved through frequency spreading implementation of 
EBMC systems compared with the EBMC implementations 
based on polyphase networks. Hence, based on the results 
of mi] and Q, we are interested in the investigation of the 
ES-EBMC in massive MIMO channels. As will be shown in 
Section HYl due to the higher frequency resolution of ES- 
EBMC, the subcarrier bands in these systems can be widened 
even more than what is reported in |I7]. Thus, the latency due 
to the synthesis and analysis filter banks can be reduced to a 
great extent compared with single antenna systems. 

In this section, we derive the mathematical foundation for 
MMSE-FSE in massive MIMO systems. In a multiuser MIMO 
system similar to the one in 0, let Nr be the number of 
receive antennas at the base station and M the number of 
users. Each user’s mobile terminal (MT) is equipped with a 
single transmit/receive antenna. All MTs are simultaneously 
communicating with the BS in time division duplexing (TDD) 
mode. The number of BS antennas is considered much larger 
than the number of users, i.e.. Nr ^ M. The users communi¬ 
cate with the BS using CMT modulation. The users’ signals 
are distinguished by their respective channel gains between 
each MT antenna and the BS antennas. The channel gain 
vectors for different users are assumed to be independent with 
respect to each other. Recalling the structure of the CMT 
receiver from Section |II] the output at the bin of the EFT 
block from the f*** user at all the receive antennas can be 
stacked into an x 1 vector rf = rfh.i{£) where is the 
channel gain vector whose elements are the gains between the 
MT antenna and the BS antennas at the frequency bin. 
The coefficient rf is the received signal at the bin of the 
output of the EFT block at the receiver from user f in an ideal 
channel. The vectors v\ are the received signals from different 


users that contribute to form the received signal at the BS. 
Accordingly, the received signal vector can be written as 

M-1 

= X! ^ 

where H* = [h^O),..., h*(M - 1)], n = [r°,..., 
and Vi is an x 1 vector that contains the EFT samples of 
the additive white Gaussian noise (AWGN) at the bin of 
the output of the EFT blocks at different receive antennas. 

In order to estimate different users’ transmitted symbols, 
a number of equalization techniques can be utilized; namely, 
matched filter (ME) and MMSE. However, as it is noted in 
Q, in realistic situations when the number of BS antennas is 
finite, ME equalization incurs some performance degradation 
due to the residual multiuser interference. In this paper, our 
discussion will be limited to the MMSE-FSE which is known 
to be a superior technique. 

In MMSE-FSE for massive MIMO systems, similar to 
single antenna ES-EBMC systems, equalization needs to be 
performed before despreading the output of the EFT blocks. 
Therefore, the MMSE linear combining aims at estimating the 
vector Yi that contains the equalized samples of different users 
at the output of their EFT blocks. The MMSE solution of dH 
is optimal as it maximizes the signal-to-interference-plus-noise 
(SINR), ll24l . Following the same line of derivations as in [T], 
the MMSE estimates of r^’s can be obtained as 

G=Wff„ (5) 

where the coefficient matrix -f 

contains the optimal MMSE filter tap weights for different 
users in its columns. It is worth mentioning that the ele¬ 
ments of the noise vector are assumed to be independent 
and identically distributed Gaussian random variables with 
variances of trj and E [v^vf ] = After the MMSE 

estimation of the vectors for i = ,N — frequency 
despreading can be separately applied for each user. To this 
end, let R = [fo,..., fjv_i]^. Consequently, the MMSE 
estimates of the transmitted data symbols of different users, 
i.e., S£ for £ = 0, ..., M—1 can be obtained through frequency 
despreading 


S = 


(6) 



where the Lx M matrix S contains the MMSE estimation of 
the transmitted data symbols of different users, i.e., s^’s, on 
its columns. 

Based on equations (|4|i to (|6]l and the assumption of having 
a flat fading channel gain over each frequency bin of DFT, 
the output SINR of the CMT receiver with MMSE channel 
equalization for the user at subcarrier i can be calculated 
as follows. 


SINR 


[i,i) __ 

MMSE ~ 


3(^.0 


3(^.0 ^ 


(7) 


where and are signal and interference-plus-noise 

powers, respectively. Pg ’ and Pj ’ can be obtained using 
equations ® and (|9]l that are shown on the top of the previous 
page where 'Wi^snijn) and s(f(TO) represent the real parts of 
the TO**' columns of the matrices and H^, respectively. 
Similarly, the subscript SJ in the equations ® and ® shows 
the imaginary part of the corresponding vectors. It is worth 
mentioning that equation O is used as a benchmark to 
investigate the channel flatness assumption in Section |IV] 


IV. Numerical Results 

In this section, the theoretical developments of the paper are 
analyzed and corroborated through numerical results. We use 
the results of m as the basis in order to evaluate the signal 
processing power of the ES method in the context of massive 
MIMO. It is worth noting that all of our simulations are based 
on a sample set of channel responses generated according 
to the SUI-4 channel model proposed by the IEEE802.16 
broadband wireless access working group, ll25l . Additionally, 
the channels between different users and different antennas are 
considered independent with respect to each other. 

In the first set of simulations (Eig. |2|i, we study the effect of 
increasing the number of antennas on the signal-to-interference 
ratio (SIR) performance. In Eig. |2] a single-user in a noise- 
free channel is considered to investigate the self-equalization 
property in ES-EBMC based massive MIMO systems. To 
evaluate the robustness of the ESE to the subcarrier width, 
we have repeated the simulations for three cases of 8, 16, and 
32 subcarriers and different number of receive antennas at 
the base station. More specifically, for the total bandwidth of 
2.8 MHz, we consider the subcarrier widths of 350, 175, and 
87.5 kHz, respectively. SIRs are evaluated at all the subcarrier 
channels. Note that in each curve, the number of points along 
the normalized frequency axis is equal to the number of 
subcarriers, L. As these results highlight, we can achieve a 
significant performance improvement even for the subcarrier 
spacings as large as 350 kHz. 

The effectiveness of these results will be more clear if 
we compare the ESE technique with the conventional PEN 
algorithm. Eig. |3] compares the SIR performance of the PPN- 
EBMC with FS-FBMC system. A single-tap equalizer per 
subcarrier is used in the PPN-FBMC structure. From Fig. |3l 
it can be understood that the SIR improvement of higher than 
30 dB can be achieved through FS-FBMC in comparison with 




Eig. 2. |(a)||(b)| and |(c)| compare the SIR performance for the ca.ses of 8,16 
and 32 subcarriers, respectively, with different number of receive antennas. 
Nr. 
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Fig. 3. SIR for the case of having L = 16, K = 4, M = 1 and different Fig. 4. SINK evaluation of the MMSE linear combining for the case of 
number of BS antennas. L = 16, K = A, and M = 6. The SNR at the receiver input is —1 dB. 


the PPN-FBMC structure. This means that higher order mod¬ 
ulation schemes can be utilized in FS-FBMC based massive 
MIMO systems compared with their PPN-based counterpart. 
The same as the previous case, the total bandwidth is fixed to 
2.8 MHz, the number of subcarriers L — 16 and overlapping 
factor of K = A are considered. This results in the subcarrier 
spacing of 2800/16 = 175 kHz. This subcarrier spacing is 
relatively wide. In fact, it is 175/15 ~ 12 times larger than 
the subcarrier spacing of the OFDM-based systems (e.g. IEEE 
802.16 and LTE). 

Next, we consider a multiuser scenario with M = 6 users, 
L = 16 subcarriers and K = 4. With the assumption of having 
perfect power control for all the users, the signal-to-noise-ratio 
(SNR) of SNRin = — 1 dB at the input of BS antennas is 
considered. The output SINR can be obtained using SNRjn + 
lOlogj^piVr in dB. Eig. |4] compares the SINR values from 
simulation for two cases of = 128 and = 64 with 
the theoretical SINR values that are calculated based on the 
assumption of having a flat fading channel per subcarrier band 
that is calculated from equation (|7]i- As the figure depicts, 
the theoretical derivations match perfectly with the simulation 
results. This emphasizes that the channel is flattened over all 
the subcarrier bands, thanks to the self-equalization property. 
The subcarrier spacing that achieves the same SINR as that 
of (|7|i is 4 times larger than what was suggested in 121 which 
brings a significant improvement in terms of shortening the 
length of the prototype filter and hence reducing the latency 
that is imposed by the transmit and receive filter banks. 

V. Conclusion 

In this paper, the idea of frequency spreading EBMC was ex¬ 
tended to the massive MIMO systems. In addition, an effective 
MMSE equalization scheme for EBMC-based massive MIMO 
systems was derived. It was shown that further subcarrier 


widening than what was suggested in Q is possible using ES- 
EBMC structure while gaining from self-equalization property. 
Based on the SIR and SINR performance results, higher 
order modulation techniques can be utilized in ES-EBMC 
based massive MIMO systems. Eurthermore, improvements 
in terms of bandwidth efficiency, robustness to frequency 
offsets, complexity, PAPR and latency compared with the 
polyphase based EBMC systems can be achieved using ES- 
EBMC structure. 
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